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Abstract 

We  present  and  illustrate  a  new  technique,  Image  Correlation  Supersampling  (ICS),  for  resampling  volume 
data  that  are  undersampled  in  one  dimension.  The  resulting  data  satisfies  the  sampling  theorem,  and, 
therefore,  many  visualization  algorithms  that  assume  the  theorem  is  satisfied  can  be  applied  to  the  data. 
Without  the  supersampling  the  visualization  algorithms  create  artifacts  due  to  aliasing. 

The  assumptions  made  in  developing  the  algorithm  are  often  satisfied  by  data  that  is  undersampled 
temporally.  Through  this  supersampling  we  can  completely  characterize  phenomena  with  measurements  at 
a  coarser  temporal  sampling  rate  than  would  otherwise  be  necessary.  This  can  save  acquisition  time  and 
storage  space,  permit  the  study  of  faster  phenomena,  and  allow  their  study  without  introducing  aliasing 
artifacts. 

The  resampling  technique  relies  on  a  priori  knowledge  of  the  measured  phenomenon,  and  applies,  in 
particular,  to  scalar  concentration  measurements  of  fluid  flow.  Because  of  the  characteristics  of  fluid  flow,  an 
image  deformation  that  takes  each  slice  image  to  the  next  can  be  used  to  calculate  intermediate  slice  images 
at  arbitrarily  fine  spacing.  We  determine  the  deformation  with  an  automatic,  multi-resolution  algorithm. 

1  Introduction 

Many  visualization  tools  assume  that  they  can  interpolate  samples  in  an  image  or  volume  to  create  a  continu¬ 
ous  function  over  a  region.  Examples  include  marching  cubes  [Lorensen  and  Cline,  1987],  which  tri-linearly 
interpolates  sample  values  to  find  isosurfaces,  and  standard  image  supersampling  [Foley  et  al.,  1990],  which 
uses  various  interpolation  techniques  to  calculate  pixel  values  between  known  ones. 

If  the  assumption  that  samples  can  be  interpolated  is  not  valid,  the  visualization  algorithms  will  produce 
results  with  artifacts.  The  assumption  is  equivalent  to  assuming  that  sampled  data  satisfies  the  sampling 
theorem  [Oppenheim  et  al.,  1983].  The  sampling  theorem  states  that  regularly  spaced  point  samples  of  a 
function  can  be  used  to  exactly  reconstruct  the  function  if  and  only  if  the  function  has  no  content  above  half 
the  sampling  frequency.  While  some  types  of  data  satisfy  the  sampling  theorem,  others  do  not.  When  a 
function  is  not  properly  sampled  the  samples  cannot,  in  general,  be  interpolated  to  reconstruct  it.  Figure  1 
shows  an  example  where  artifacts  are  introduced  by  interpolative  supersampling  of  an  undersampled  image. 
High  frequency  information  encoded  in  the  sample  values  aliases  into  lower-frequency  artifacts  in  the 
supersampled  data. 

1.1  Creating  Properly  Sampled  Data:  Image  Correlation  Supersamping 

We  have  developed  a  technique  for  converting  certain  sampled  data  that  do  not  satisfy  the  sampling  theorem 
into  sampled  data  that  do.  Our  primary  application  involves  temporally  supersampling  time- varying  images 
of  flowing  fluid  where  the  temporal  sampling  rate  is  too  low,  but  the  technique  may  also  be  applied  to  volume 
data  where:  1)  one  dimension  is  undersampled,  2)  a  parameterized  deformation  characterizes  the  change 
from  one  slice  image  to  the  next,  and  3)  the  parameters  of  the  deformation  between  two  slices  can  be  found. 
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Figure  1 :  An  example  of  the  impact  of  interpolative  methods  for  supersampling  undersampled  data,  (a) 
A  properly  sampled  image  of  a  shape  with  a  Gaussian  cross  section,  (b)  shows  supersampling  of  an 
undersampled  version  of  the  same  object.  The  undersampling  and  supersampling  are  in  the  vertical  direction. 
Note  the  bumpy  shapes  introduced. 


The  algorithm  we  apply,  which  is  described  in  Section  2  and  in  more  detail  in  [Gornowicz,  1997],  uses  a 
multi-resolution  deformation  to  characterize  the  change  from  one  slice  image  to  the  next  and  automatically 
optimizes  the  deformation  parameters  at  each  level  of  resolution  to  find  the  appropriate  deformation.  Once 
the  deformation  has  been  found,  it  is  used  to  calculate  intermediate  images  at  a  sufficient  sampling  rate  to 
satisfy  the  sampling  theorem. 

1.2  Related  Work 

This  work  has  been  motivated  in  part  by  metamorphosis,  or  morphing,  techniques  for  both  images 
[Beier  and  Neely,  1992]  and  volumes  [Lerios  et  al.,  1995].  Both  of  these  techniques  use  carefully  crafted 
user  input  to  define  a  deformation  that  takes  one  image  or  volume  into  another.  For  intermediate  frames 
they  both  dissolve  from  one  partially  deformed  image  or  volume  into  the  second.  Our  technique  differs 
primarily  in  that  the  deformation  we  calculate  is  physically  motivated  by  the  equations  of  fluid  flow,  and  is 
found  automatically  rather  than  being  interactively  specified. 

A  similar  deformation-based  technique  [Ruprecht  and  Muller,  1994]  defines  correspondence  points  on 
different  magnetic  resonance  images  and  calculates  a  deformation  by  attempting  to  deform  the  images  to 
satisfy  the  correspondences.  Our  deformation  has  a  more  physical  basis,  and  is  parameterized  and  discovered 
differently. 

Video  compression  techniques  solve  a  similar  problem.  By  using  a  predictive  model  of  the  deformation 
of  one  video  frame  to  the  next  they  can  encode  only  the  differences  between  a  deformed  version  of  an 
earlier  frame  and  the  current  frame  [Huang  et  al.,  1996].  If  the  deformation  is  accurate,  then  the  differences 
will  be  small  and  the  resulting  compression  high.  Because  a  failure  in  the  deformation  leads  only  to  less 
effective  compression,  these  techniques  do  not  fail  at  representing  the  information  if  they  are  not  robust, 
their  compression  rates  degrade.  Although  our  technique  is  limited  to  certain  types  of  image  data,  it  robustly 
supersamples  them  to  satisfy  the  sampling  theorem,  while  the  video  compression  techniques  would  not. 

This  work  is  based  on  Image  Correlation  Velocimetry  (ICV)  [Tokumaru  and  Dimotakis,  1995],  an  auto¬ 
matic  technique  for  extracting  velocity  information  from  a  series  of  images  of  concentration  data.  That  orig¬ 
inal  work  was  subsequently  extended  with  a  multi-resolution  basis  for  the  deformation  [Gornowicz,  1997]. 
The  velocity  information  serves  as  the  deformation  in  the  two-step  metamorphosis  process. 
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Other  related  techniques  for  the  computational  analysis  of  motion  from  a  sequence  of  images  include  par¬ 
ticle  image  velocimetry  [Adrian,  1991]  [Willert  and  Gharib,  1991],  variational  approaches  [Zhou  et  al.,  1995] 
[Szeliski  and  Shum,  1996],  and  differential  techniques  [Horn  and  Schunck,  1981]  [Willick  and  Yang,  1991] 
[Dahm  et  al.,  1992]).  A  review  of  so-called  optical  flow  techniques  is  given  by  [Barron  et  al.,  1994]. 

1.3  Roadmap 

Section  2  describes  the  algorithm  and  is  followed  by  a  description  of  the  results  produced  by  applying  the 
algorithm  to  both  simulated  and  acquired  data.  A  discussion  of  the  algorithm  and  results  is  in  Section  5, 
with  conclusions  immediately  after. 

2  Image  Correlation  Supersampling 

The  Image  Correlation  Supersampling  (ICS)  algorithm  consists  of  two  parts,  as  do  many  other  morphing 
techniques.  In  the  first  part  a  deformation  or  displacement  field  is  defined  that  transforms  one  image  into  the 
next.  In  the  second  part,  a  weighted  sum  of  the  partially  displaced  images  is  used  to  generate  intermediate 
images  at  a  sampling  rate  sufficient  to  satisfy  the  sampling  theorem.  This  rate  can  be  determined  from 
knowledge  of  the  phenomena  being  measured  or  by  estimating  the  rate  conservatively  and  examining  the 
resulting  samples  to  calculate  a  more  accurate  and  lower  rate.  For  fluid  experiments,  a  maximum  velocity 
is  frequently  available  and  can  be  used  to  bound  the  maximum  sampling  rate. 

In  this  section  we  first  present  the  physically-based  motivation  for  ICS,  and  then  discuss  the  two  parts 
of  the  algorithm. 

2.1  Physical  Motivation  for  Deformation 

For  2D  flows  where  a)  diffusive  effects  are  negligible  (e.g.  high  Schmidt  number  or  short  time  periods),  and 
b)  the  concentration  is  proportional  to  image  intensity,  the  scalar  transport  equation  applies  to  the  intensity 
image  7(x,  t),  and 


where  u  is  the  velocity  [Batchelor,  1967]. 
displacement  field  £  =  ru ,  then 


d  d 


(i) 


With  r  a  small  step  in  time,  and  with  the  definition  of  the 


4,+(4/^°  <2) 

For  two  concentration  images,  7(x,  t)  and  7(x,  t  +  r),  7(x,  t  +  r)  can  be  written  as  a  Taylor  expansion  around 
x  and  t. 


I(x  +  Z,t+T)  =  I(x.  t)  +  t  t)  +  £  ■  7— /(x.  t)  +  H.O.T.  (3) 

For  small  £  and  r,  the  higher  order  terms  can  be  neglected.  The  sum  of  the  first  order  terms  in  £  and  r  are 
also  negligible  from  Equation  2. 

2.2  Image  Deformation 

ICS  seeks  the  displacement  field,  ^(x),  that  transforms  one  images  into  the  next.  We  define  7i(x)  =  7(x,  t  +  r) 
and  7q(x)  =  7(x,  t). 
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I[(x  +  £(x))  ~  70(x) 


(4) 


<£(x)  is  a  2-vector  field  over  the  spatial  coordinates,  x,  where  our  image  is  defined,  <^(x)  represents  a 
displacement,  or  deformation,  of  the  spatial  coordinates,  and  hence  of  the  image.  Given  the  two  scalar  field 
images  7q(x)  at  time  t  and  7i(x)  at  time  t  +  r,  ICS  seeks  a  displacement  field  ^(x)  that  minimizes  the  cost 
function  /(£)  given  by  the  square  of  the  differences  in  intensity  integrated  over  the  correlation  domain,  fl. 

min  m  =  min  f  +  ^  _  /q(jc)]2  (5) 

£  £  Jn 

The  displacement  field,  ^(x),  is  implemented  using  a  multi-resolution  hierarchy  of  tri-cubic  B-splines,  and 
the  minimization  of  /(£)  is  implemented  with  a  multi-dimensional  gradient  approach.  More  details  are 
available  in  [Gornowicz,  1997]. 

2.3  Supersampling 

Intermediate  slice  images,  subscripted  by  0  <  a  <  1,  w  are  computed  from  the  surrounding  images  and  the 
displacement  field,  ^(x),  using 


7a(x)  =  (1  -  a)70(x  -  aO  +  a7i(x  +  (1  -  a)£)  (6) 

where  a  varies  from  0  to  1  as  the  image  transitions  from  7q  to  I\ .  The  equation  linearly  blends  image  7q  to 
7i,  with  only  image  7o  contributing  when  a  =  0.0  and  only  image  I\  contributing  when  a  =  1.0  Image  7o  is 
linearly  displaced;  it  is  undeformed  when  a  =  0.0  and  fully  displaced  when  a  =  1.0.  I\  is  linearly  inversely 
displaced;  it  is  fully  displaced  when  a  =  0.0  and  undeformed  when  a  =  1 .0. 

3  Imaging 

In  this  section  we  describe  the  two  volume  images  created  and  the  one  acquired  to  test  and  illustrate  our 
supersampling  technique.  In  the  two  simulated  volumes  the  sampling  theorem  is  satisfied  in  two  of  the 
three  directions.  In  the  third  direction  the  sampling  theorem  is  not  satisfied-the  sampling  kernel,  or  point 
spread  function,  is  too  narrow;  equivalently,  the  samples  are  spaced  too  far  apart.  However,  the  behavior  of 
the  simulated  phenomena  in  the  regions  between  slices  is  well  characterized  as  a  deformation  of  the  nearby 
slices,  and  thus  satisfies  our  assumptions.  In  the  acquired  data  the  two  spatial  sampling  rates  are  again 
sufficient  to  satisfy  the  sampling  theorem,  but  the  temporal  sampling  rate  is  not. 

3.1  Translated  Disk 

For  the  first  simulated  dataset  a  2D  normal  distribution  translates  from  the  lower  left  toward  the  upper  right 
of  the  image  (see  Figure  2).  The  dataset  consists  only  of  the  two  images,  and  the  deformation  is  a  simple 
translation. 

3.2  Translated,  Rotated,  Strained  Images 

For  the  second  simulated  dataset  an  image  of  two  regions  of  band-limited  noise  are  translated,  rotated,  and 
strained  in  two  dimensions.  Figure  3  shows  three  frames  from  the  dataset.  These  deformations  are  exact 
solutions  to  the  Navier-Stokes  equations  for  2D  flow. 
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Figure  2:  Frames  from  simulated  dataset  of  translating  disk. 


Figure  3:  Frames  from  a  simulated  2D  flow  experiment. 


Figure  4:  Frames  from  measurements  of  turbulent  fluid  flow.  Time  increase  for  images  toward  the  right. 
The  overall  flow  within  each  image  is  toward  the  left. 


3.3  Turbulent  Jet  Data 

The  experimentally  measured  data,  shown  in  Figure  4,  consists  of  images  of  a  liquid-phase,  axisymmetric, 
turbulent  jet,  with  Re  -  2250,  discharging  into  a  large  quiescent  reservoir  of  square  cross-section.  The  flow 
was  established  by  pressurizing  the  jet  plenum  in  a  blow-down,  nearly-constant-pressure  configuration  from 
two  gas  cylinders.  The  jet  nozzle  had  an  exit  diameter  of  2.54  mm,  and  the  discharge  tank  had  internal 
dimensions  of  0.85  x  0.85  x  1.6m  (335  x  335  x  625  diameters).  Rectangular  glass  windows  on  all  four 
sides  of  the  tank,  in  combination  with  a  10  inch  diameter  glass  window  at  the  bottom  of  the  tank,  provide 
optical  access  to  the  flow  facility. 

A  laser  sheet  containing  the  axis  of  the  jet  was  formed  by  passing  the  beam  from  a  frequency-doubled 
Nd:YAG  laser  pulsed  at  10  Hz  through  a  —6.4mm  focal  length  cylindrical  lens.  A  long  focal  length 
(f  -  +1.0m)  cylindrical  lens  oriented  90  degrees  with  respect  to  the  first  lens  ensured  that  the  laser  sheet 
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Figure  5:  2D  slices  of  simulated  3D  data,  with  the  undersampled  axis  vertical  and  one  of  the  sufficiently 
sampled  axes  horizontal,  (a)  shows  Image  Correlation  Supersampling  (ICS),  and  (b)  interpolative  super¬ 
sampling.  In  the  3D  data  a  bright  disk  translates  from  left  to  right,  as  is  reconstructed  correctly  by  ICS  in 
(a).  The  disk  is  incorrectly  split  into  two  features  in  (b). 


(a)  (b) 


Figure  6:  (a)  Image  Correlation  Supersampling  (ICS),  and  (b)  interpolative  supersampling.  The  images 
show  an  isosurface  calculated  from  the  3D  data.  In  the  3D  data  a  bright  disk  translates  from  left  to  right,  as 
is  reconstructed  correctly  by  ICS  in  (a).  The  disk  is  incorrectly  split  into  two  features  in  (b). 


reached  its  waist,  and  remained  thin,  in  the  region  of  interest.  The  532  nm  line  from  the  laser  was  used 
to  excite  a  fluorescent  dye  that  was  pre-mixed  in  weak  concentration  with  the  jet  plenum  fluid.  The  laser- 
induced  fluorescence  field  was  recorded  at  a  rate  of  10  frames/second  on  a  12-bit  camera  that  utilizes  a 
1024  x  1024  pixel  Cassini  CCD.  A  sequence  of  42  images  of  a  125  X  125  diameter  region  centered  on  the 
jet  axis  at  a  downstream  distance  of  225  diameters  from  the  jet  nozzle  was  recorded.  For  this  experiment, 
the  in-plane  spatial  resolution  and  the  out-of-plane  spatial  (laser-sheet  thickness)  resolution  was  sufficient 
to  resolve  the  diffusion  scale. 
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Figure  7:  2D  slices  of  3D  data  measuring  turbulent  flow,  with  the  undersampled  temporal  axis  vertical 
and  one  sufficiently  sampled  spatial  axis  horizontal,  (a)  Image  Correlation  Supersampling  (ICS),  and  (b) 
interpolative  supersampling,  (b)  shows  significant  aliasing  introduced  by  the  supersampling.  The  aliasing 
is  not  apparent  in  (a),  (c)  is  a  gold  standard  for  the  supersampled  data. 


(a)  (b)  (c) 


Figure  8:  An  isosurface  calculated  from  simulated  2D  fluid  flow,  (a)  Image  Correlation  Supersampling 
(ICS),  and  (b)  interpolative  supersampling,  (c)  shows  an  isosurface  calculated  from  data  correctly  sampled 
at  high  resolution. 


4  Results 

The  image  correlation  supersampling  technique  has  been  applied  to  the  datasets  described  above.  We 
show  the  results  of  some  common  visualization  techniques  applied  to  the  data  both  with  and  without  ICS 
supersampling. 

The  first  example  of  Image  Correlation  Supersampling  is  shown  in  Figures  5  and  6.  The  original  2D 
image  data  for  this  example  represent  a  translating  disk  shape  and  are  shown  in  Figure  2.  Figure  5  shows  2D 
slices  taken  perpendicular  to  the  original  images,  with  one  spatial  dimension  on  the  horizontal  axis  and  time 
on  the  vertical  axis.  Supersampling  was  accomplished  via  ICS  in  5(a)  and  via  interpolative  supersampling 
in  5(b).  Note  that  ICS  supersampling  correctly  maintains  the  integrity  of  the  feature  translating  through  the 
3D  volume,  whereas  the  feature  is  split  into  two  parts  in  (b).  Figure  6  shows  an  isosurface  created  from  the 
same  data.  In  6(a),  the  correct  translation  of  the  disk  leads  to  a  rod,  as  expected,  while  in  6(b)  the  rod  is 
incorrectly  divided  into  two  separate  features. 

The  second  example  of  ICS  is  shown  in  Figures  7  and  8.  The  original  2D  image  data  for  this  example 
represent  two  deforming  regions  of  band-limited  noise  (see  Figure  3).  Figure  7  shows  2D  slices  taken 
perpendicular  to  the  original  images,  with  one  spatial  dimension  on  the  horizontal  axis  and  time  on  the 
vertical  axis.  Supersampling  was  accomplished  via  ICS  in  7(a)  and  via  interpolative  supersampling  in  7(b). 
Note  that  ICS  supersampling  correctly  maintains  the  integrity  of  features  translating  through  the  3D  volume, 
whereas  the  features  are  distorted  in  7(b).  Figure  7(c)  shows  the  data  directly  sampled  at  high  temporal 
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(a)  (b) 


Figure  9:  An  isosurface  calculated  from  3D  imaging  data  of  turbulent  fluid  flow,  (a)  Image  Correlation 
Supersampling  (ICS),  and  (b)  interpolative  supersampling. 


Figure  10:  (a)  Image  Correlation  Supersampling  (ICS),  and  (b)  interpolative  supersampling.  The  images 
show  2D  slices  of  3D  data  measuring  turbulent  flow,  (b)  shows  significant  aliasing  introduced  by  the 
supersampling.  The  aliasing  is  not  apparent  in  (a). 


resolution  from  the  analytic  expression  for  the  fluid  flow.  It  is  a  gold  standard  for  the  supersampled  data. 
Figure  8  shows  an  isosurface  created  from  the  same  data.  In  8(a),  the  correct  deformation  of  the  noise 
images  leads  to  two  slightly  twisted  tubes,  as  expected,  while  in  8(b)  the  tubes  are  split  into  pieces  and  are 
significantly  distorted.  8(c),  a  gold  standard  for  this  test  case,  shows  what  the  isosurface  should  look  like 
given  perfect  sampling. 

We  have  also  applied  the  ICS  algorithm  to  data  imaged  from  turbulent  flow.  The  data  for  this  example  are 
shown  in  Figure  4,  and  are  sampled  at  10  frames/second.  The  results  are  supersampled  at  200  frames/second, 
a  rate  determined  to  be  sufficient  from  the  velocity  of  the  jet.  The  45°  angle  of  features  in  Figure  10  further 
suggests  that  the  rate  is  sufficient  because  the  temporal  sampling  of  the  image  is  similar  to  the  spatial 
sampling,  where  the  sampling  rate  is  known  to  be  sufficient. 

Figure  9(a)  shows  an  isosurface  of  concentration  data  measured  in  a  plane.  Time  is  stacked  in  the  third 
dimension.  The  jet  flows  towards  the  left,  and  time  moves  from  top  to  bottom  in  the  box.  Note  that  structures 
moving  through  time  trace  out  smooth,  slanted  surfaces,  with  steep  surfaces  representing  slow  structures, 
and  shallow  surfaces  fast  structures.  In  Figure  9(b)  interpolative  supersampling  introduces  many  artifacts, 
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analogous  to  those  shown  in  earlier  figures,  which  are  incorrect  and  which  mislead  the  viewer  about  the 
behavior  of  the  fluid.  Note  the  ridges  on  the  left-most  surface  in  the  image,  and  the  three  separate  shapes 
at  the  far  right  corresponding  to  a  single  rod  in  Figure  9(a).  Figure  10  shows  a  2D  section  of  the  3D  data. 
Note  the  artifacts  in  10(b)  created  by  aliasing.  They  are  absent  in  10(a),  the  ICS  image. 

5  Discussion 

We  have  made  some  observations  about  ICS  and  its  application  while  working  to  apply  it  to  the  example 
datasets. 

5.1  Interpreting  Image  Displacements  as  Velocity 

The  ICV  algorithm  [Tokumaru  and  Dimotakis,  1995]  interprets  the  displacement  field  that  it  calculates  as 
a  velocity.  For  some  cases,  in  particular  where  image  intensities  are  constant  within  a  region,  the  velocity 
is  unconstrained,  and  so  there  are  many  possible  solutions.  This  ambiguity  leads  to  uncertainty  in  the 
velocities.  For  ICS,  however,  the  ambiguity  in  the  displacement  is  canceled  out  in  the  supersampling  step 
because  the  displacement  ambiguity  lies  completely  within  regions  of  constant  intensity. 

5.2  Out-of-plane  Motion 

The  physical  basis  for  using  a  displacement  field  to  represent  fluid  behavior  between  images  is  based  on  2D 
flow  within  the  images.  In  cases  where  there  is  flow  out  of  the  imaging  plane  the  changes  are  not  correctly 
accounted  for  in  intermediate  images.  Some  artifacts,  then,  remain  in  the  supersampled  data. 

The  jet  data  that  we  measured  (Figure  4)  does  have  fluid  flow  components  out  of  the  imaging  plane. 
Because  of  the  orientation  of  the  jet  along  the  imaging  plane,  the  out-of-plane  velocity  components  are 
smaller  than  ^  the  in-plane  components,  and  so  are  not  as  important.  There  are  also  significant  portions 
of  the  dataset  where  the  out-of-plane  velocities  are  negligible.  Even  with  this  out-of-plane  velocity,  ICS 
reduces  the  artifacts  due  to  high-frequency  in-plane  velocities,  and  does  not  increase  the  artifacts  in  regions 
of  out-of-plane  motion.  We  feel  that  the  costs,  then,  are  minimal  and  the  benefits  significant,  although  any 
artifacts  that  remain  must  be  carefully  interpreted. 

To  address  this  problem  we  propose  two  possibilities.  The  first  is  higher-speed  acquisitions,  although  the 
associated  cost  is  high.  The  second  is  collecting  volume  data,  so  that  a  displacement  field  can  be  calculated 
that  takes  into  account  the  3D  flow. 

5.3  Time-sampled  Volume  Data:  the  Fourth  Dimension 

A  logical  extension  of  this  algorithm  would  apply  to  4D  data-volumes  that  are  sampled  over  time.  The 
deformation  or  displacement,  £,  would  be  extended  to  3D.  This  extension  would  slow  the  optimization 
process,  but  would  be  particularly  useful  in  this  context  because  it  would  help  reduce  the  temporal  sampling 
rate  and  hence  reduce  the  long  acquisition  times  and  large  data  size  for  these  datasets. 

5.4  Other  Application  Areas 

With  an  extension  to  volumes,  and  an  appropriate  modification  to  the  type  of  displacement  of  deformation 
allowed,  this  work  should  apply  to  several  medical  imaging  areas.  MRI  time- series  are  often  undersampled 
temporally,  and  so  do  not  sample  the  motion  of  the  subject  well.  The  body,  however,  is  generally  displacing 
from  volume  image  to  volume  image,  and  the  displacement  should  be  calculatable  and  usable  to  supersample 
the  data  temporally. 
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Similarly,  developmental  studies  of  animals  track  the  changes  in  shape  of  different  features.  If  these  can 
be  modeled  as  displacements  within  volume  images,  with  expansion,  contraction,  and  deformation  of  the 
shapes  taken  into  account,  then  a  relatively  small  number  of  volume  images  of  a  developing  animal  might 
be  used  to  reconstruct  a  continuous  representation  of  the  animal  during  development. 

6  Conclusions 

The  technique  that  we  have  presented,  and  others  with  a  similar  physical  motivation,  can  be  used  to  reduce 
the  high  sampling  rates  required  of  acquired  data  to  satisfy  the  sampling  theorem  by  post-processing  the 
undersampled  data.  Computation  and  visualization  tools  that  assume  the  sampling  theorem  is  satisfied 
can  then  be  applied  to  the  measurements.  The  savings  in  collection  time,  higher- speed  hardware,  and 
storage  can  be  significant.  The  acquired  dataset  illustrating  the  algorithm  would  have  had  to  be  sampled  20 
times  as  fast  to  achieve  similar  results  without  ICS.  To  collect  at  that  rate  would  have  required  significant 
hardware  upgrades  for  both  the  excitation  laser  and  for  the  camera  and  data  acquisition  hardware-upgrades 
that  approach  and  may  even  exceed  the  capabilities  of  available  technology.  With  ICS,  existing  lower-cost 
hardware  was  sufficient,  and  the  behavior  of  the  jet  can  still  be  studied  with  standard  visualization  tools. 
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